Turbulent Coronal Heating Mechanisms: Coupling of Dynamics and Thermodynamics

Context. Photospheric motions shuffle the footpoints of the strong axial magnetic field that threads coronal loops giving rise to turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets where energy is deposited at small-scales and the heating occurs. Previous studies show that current sheets thickness is orders of magnitude smaller than current state of the art observational resolution (~700 km). Aim. In order to understand coronal heating and interpret correctly observations it is crucial to study the thermodynamics of such a system where energy is deposited at unresolved small-scales. Methods. Fully compressible three-dimensional magnetohydrodynamic simulations are carried out to understand the thermodynamics of coronal heating in the magnetically confined solar corona. Results. We show that temperature is highly structured at scales below observational resolution and nonhomogeneously distributed so that only a fraction of the coronal mass and volume gets heated at each time. Conclusions. This is a multi-thermal system where hotter and cooler plasma strands are found one next to the other also at sub-resolution scales and exhibit a temporal dynamics.Comment: A&A Letter, in pres

Turbulence, Energy Transfers and Reconnection in Compressible Coronal Heating Field-line Tangling Models

MHD turbulence has long been proposed as a mechanism for the heating of coronal loops in the framework of the Parker scenario for coronal heating. So far most of the studies have focused on its dynamical properties without considering its thermodynamical and radiative features, because of the very demanding computational requirements. In this paper we extend this previous research to the compressible regime, including an energy equation, by using HYPERION, a new parallelized, viscoresistive, three-dimensional compressible MHD code. HYPERION employs a Fourier collocation -- finite difference spatial discretization, and uses a third-order Runge-Kutta temporal discretization. We show that the implementation of a thermal conduction parallel to the DC magnetic field induces a radiative emission concentrated at the boundaries, with properties similar to the chromosphere--transition region--corona system.Comment: 4 pages, 4 figures, Solar Wind 12 proceedings (in press

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow

Modeling the Galactic Center Nonthermal Filaments as Magnetized Wake

We simulate the Galactic Center nonthermal filaments as magnetized wakes formed dynamically from amplification of a weak (tens of $\mu$G) global magnetic field through the interaction of molecular clouds with a Galactic Center wind. One of the key issues in this cometary model is the stability of the filament against dynamical disruption. Here we show 2-dimensional MHD simulations for interstellar conditions that are appropriate for the Galactic Center. The structures eventually disrupt through a shear driven nonlinear instability but maintain coherence for lengths up to 100 times their width as observed. The final instability, which destroys the filament through shredding and plasmoid formation, grows quickly in space (and time) and leads to an abrupt end to the structure, in accord with observations. As a by-product, the simulation shows that emission should peak well downstream from the cloud-wind interaction site.Comment: postscript file, 7 figs (included); Accepted for publication in ApJ (Part 1

Nonlinear Dynamics of the Parker Scenario for Coronal Heating

The Parker or field line tangling model of coronal heating is studied comprehensively via long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy. In physical space this corresponds to a magnetic topology where magnetic field lines are barely entangled, nevertheless current sheets (corresponding to the original tangential discontinuities hypothesized by Parker) are continuously formed and dissipated. Current sheets are the result of the nonlinear cascade that transfers energy from the scale of convective motions ($\sim 1,000 km$) down to the dissipative scales, where it is finally converted to heat and/or particle acceleration. Current sheets constitute the dissipative structure of the system, and the associated magnetic reconnection gives rise to impulsive ``bursty'' heating events at the small scales. This picture is consistent with the slender loops observed by state-of-the-art (E)UV and X-ray imagers which, although apparently quiescent, shine bright in these wavelengths with little evidence of entangled features. The different regimes of weak and strong MHD turbulence that develop, and their influence on coronal heating scalings, are shown to depend on the loop parameters, and this dependence is quantitatively characterized: weak turbulence regimes and steeper spectra occur in {\it stronger loop fields} and lead to {\it larger heating rates} than in weak field regions.Comment: 22 pages, 18 figures, uses emulateapj, for mpeg file associated to Figure 17e see (temporarily) http://www.df.unipi.it/~rappazzo/arxiv/jfl.mpg, ApJ, in pres

The effect of the solar corona on the attenuation of small-amplitude prominence oscillations. I. Longitudinal magnetic field

Context. One of the typical features shown by observations of solar prominence oscillations is that they are damped in time and that the values of the damping times are usually between one and three times the corresponding oscillatory period. However, the mechanism responsible for the attenuation is still not well-known. Aims. Thermal conduction, optically thin or thick radiation and heating are taken into account in the energy equation, and their role on the attenuation of prominence oscillations is evaluated. Methods. The dispersion relation for linear non-adiabatic magnetoacoustic waves is derived considering an equilibrium made of a prominence plasma slab embedded in an unbounded corona. The magnetic field is orientated along the direction parallel to the slab axis and has the same strength in all regions. By solving the dispersion relation for a fixed wavenumber, a complex oscillatory frequency is obtained, and the period and the damping time are computed. Results. The effect of conduction and radiation losses is different for each magnetoacoustic mode and depends on the wavenumber. In the observed range of wavelengths the internal slow mode is attenuated by radiation from the prominence plasma, the fast mode by the combination of prominence radiation and coronal conduction and the external slow mode by coronal conduction. The consideration of the external corona is of paramount importance in the case of the fast and external slow modes, whereas it does not affect the internal slow modes at all. Conclusions. Non-adiabatic effects are efficient damping mechanisms for magnetoacoustic modes, and the values of the obtained damping times are compatible with those observed.Comment: Accepted in A&

Magnetic Effects at the Edge of the Solar System: MHD Instabilities, the de Laval nozzle Effect and an Extended Jet

To model the interaction between the solar wind and the interstellar wind, magnetic fields must be included. Recently Opher et al. 2003 found that, by including the solar magnetic field in a 3D high resolution simulation using the University of Michigan BATS-R-US code, a jet-sheet structure forms beyond the solar wind Termination Shock. Here we present an even higher resolution three-dimensional case where the jet extends for $150AU$ beyond the Termination Shock. We discuss the formation of the jet due to a de Laval nozzle effect and it's su bsequent large period oscillation due to magnetohydrodynamic instabilities. To verify the source of the instability, we also perform a simplified two dimensional-geometry magnetohydrodynamic calculation of a plane fluid jet embedded in a neutral sheet with the profiles taken from our 3D simulation. We find remarkable agreement with the full three-dimensional evolution. We compare both simulations and the temporal evolution of the jet showing that the sinuous mode is the dominant mode that develops into a velocity-shear-instability with a growth rate of $5 \times 10^{-9} sec^{-1}=0.027 years^{-1}$. As a result, the outer edge of the heliosphere presents remarkable dynamics, such as turbulent flows caused by the motion of the jet. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how this outer boundary behaves. Already, however, we can say that the magnetic field effects are a major player in this region changing our previous notion of how the solar system ends.Comment: 24 pages, 13 figures, accepted for publication in Astrophysical Journal (2004